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Abstract 

The observed pattern of neutrino mixing may be the result of a 2-3 ( fi — r) symmetry 
in the leptonic sector. We consider a two Higgs doublet model with a 2-3 symmetry 
in the down type quark and the charged lepton sector. The breaking of the 2-3 
symmetry by the strange quark mass and the muon mass leads to FCNC in the 
quark sector and the charged lepton sector that are suppressed by ^ and ^ in 
addition to the mass of the heavy Higgs boson of the second Higgs doublet. A Higgs 
boson mass of mn ~ 600 — 900 GeV can explain the deviation from standard model 
reported in several rare B decays. Predictions for other B decays are made and a 
new CP phase is predicted in Bg — Bg mixing. The lepton flavour violating decays 
r —>■ fJ'l{q)l{q) are below the experimental limits. The breaking of 2-3 symmetry in 
the lepton sector can lead to deviations of the atmospheric neutrino mixing angle 
from the maximal value by ~ 2 degrees. 
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1 Introduction 



The discovery of neutrino masses and mixing have led to many speculations about 
extension of the standard model(SM). In contrast to the mixing in the quark sector 
which is hierarchial, the neutrino mixing is large. A useful framework to under- 
stand the large neutrino mixing are models that have a leptonic /i — r interchange 
symmetry ^ 121. Several ongoing and future neutrino experiments are expected to 
provide us with insights into the physics behind neutrino masses and mixing. 

On the other hand there is lot of data now available on CP violation in the B 
system. The goal of the B factories is to test the CKM picture of CP violation and 
look for evidence of new physics. 

For some time now the B factories have been reporting several experimental hints 
of new physics. First, within the SM, the measurement of the CP phase sin 2/5 in 
B^{t) —>■ J/ipKs should be approximately equal to that in decays dominated by the 
quark-level penguin transition b — > sqq (g = u, d, s) like B^{t) (pKg, B^{t) — >• rj'Kg, 
B^[t) TT^K^, etc. However, there is a difference between the measurements of 
sin 2/3 in the b ^ s penguin dominated modes (sin 2/3 = 0.50 ± 0.06) and that in 
Bdit) J/ipKsism2p = 0.685 ± 0.032) 0111 El- Note that the sin2/? number for 
the b ^ s penguin dominated modes is the average of several modes. The effect of 
new physics can be different for different final hadronic states and so the individual 
sin 2/? measurements for the different modes are important. Second, the latest data 
on B ^ ttK decays (branching ratios and various CP asymmetries) appear to be 
inconsistent with the SM [HI El- ^ Third, within the SM, one expects no triple- 
product asymmetries in B ^ (pK* jU], but BaBar has measured such an effect at 
1.7a level ^01 • There are also polarization anomalies where the decays B^ (pK* 
and p^K* appear to have large transverse polarization amplitudes in conflict 

with naive SM expectations [51 ITU IT^. 

While these deviations certainly do not unambiguously signal new physics(NP), 
they give reason to speculate about NP explanations of the experimental data. Fur- 
thermore, it is far more compelling to find NP scenarios that provide a single solution 
to all the deviations than to look for solutions to individual discrepancies. Taking 
all these deviations seriously one is lead to certain structures of NP operators that 
can explain the present data^^. The question then is what kind of NP models 
can generate these specific operator structures at low energies. One can be more 
ambitious and ask if such models may have any connection with the physics be- 
hind neutrino mixing. In fact, the large z/^ — Ut mixing can lead to new effects in 
6 — >■ s transitions through squark mixing in certain Supersymmetric Grand Unified 
Theories P^j. However, the dominant contributions in this scenario are new QCD 
penguin amplitudes to B decays involving b ^ s transitions and such contributions 
appear to be ruled out by the present B — > Kit data[ni[7j. 

cleaner test of the SM could be provided by looking at the quasi-exclusive decays B KX 
|S] rather than the exclusive B Kir decays. 
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In this paper we present a simple two Higgs doublet model with a 2-3 interchange 
symmetry in the down quark sector like the H — t interchange symmetry in the 
leptonic sector. This model can provide a single solution to the deviations from SM 
observed in B decays and has interesting implications for the leptonic sector. It is 
believed that one of the phenomenological problems of the 2-3 symmetry model is 
the predictions = rrir or rus = rrtb. However, we will assume the 2-3 symmetry in 
the gauge basis where the mass matrix has off diagonal terms and is 2-3 symmetric. 
Diagonalizing the mass matrix will split the masses of s and b or fi and r. In fact 
we will consider the "enhanced" 2-3 symmetry where the matrix element of the 
mass matrix are invariant under any 2-3 interchange. In other words we assume 
< 2|M|2 >=< 2|M|3 >=< 3|M|2 >=< 3|M|3 >. Diagonalizing the mass matrix 
then leads to vanishing ms(m^). 

The breaking of the 2-3 symmetry is then introduced though the strange quark 
mass in the quark sector and the muon mass in the leptonic sector. The breaking 
of the 2-3 symmetry leads to flavor changing neutral currents (FCNC) in the quark 
sector and the charged lepton sector that are suppressed by ^ and ^ in addition 
to the mass of the Higgs boson of the second Higgs doublet. Additional FCNC 
effects of similar size can be generated from the breaking of the s — b symmetry in 
the Yukawa coupling of the second Higgs doublet. 

In this paper we will be interested in the simplest two Higgs doublet model with a 
2-3 symmetric Yukawa coupling that can explain the hints of new physics observed 
in several rare B decays. We, therefore, make some simplifying assumptions to 
make our model as predictive as possible. First, we assume a discrete symmetry 
involving the down quark to prevent FCNC effects in s — *• (i and b ^ d transitions. 
This allows us to satisfy constraints from measurements in the kaon system and Bd 
mixing. Second, we assume a simple ansatz for the Yukawa coupling of the second 
Higgs doublet where the Yukawa matrix is described by two real parameters and 
a universal weak phase. Finally, we assume there is no mixing among the neutral 
Higgs bosons in the model and we only consider FCNC effects generated by the 
lightest Higgs bosons of the second Higgs doublet. Possible FCNC effects from the 
other heavier neutral Higgs boson are neglected by assuming the mass of the Higgs 
boson to be sufficiently high. 

The low energy effective Hamiltonian generated by this model can explain the 
general features of the deviations from the SM seen in the various rare B decays. A 
new physics fit to the B Kit data j7] allows us to fix the parameters of the model 
including the lightest Higgs mass of the second doublet. To extract the parameters 
of the model we have to calculate the ratio of the tree amplitude in the standard 
model relative to a new physics amplitude. We use factorization to calculate the tree 
amplitude and the new physics amplitudes in the Ktt system. Since non factorizable 
effects in the tree and the new physics amplitudes are expected to be small [T31 HH] 
the use of factorization to extract the parameters of the model is reasonable. A more 
precise determination of the parameters of the model using the calculation of the 
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non leptonic amplitudes in the framework of QCD factorization will be carried 
out in a future work.l7j. 

Having fixed the parameters of the model we make predictions for the rare B 
decays B —>■ pK* ,(f)K* ,(f)Ks and rj'Kg. Here again we use factorization to calculate 
the relevant non leptonic amplitudes. The inclusion of non factorizable effects will 
change the predicted values of quantities such as polarization fractions, CP violation 
etc, but we do not expect the general pattern of new physics effects of the model 
in the various B decays to change. This is because our predictions are largely 
dependent on heavy quark theory considerations. More precise predictions for the 
various B decays in this model, including non factorizable effects, will be carried out 
in a future work[17j. 

We also make predictions for B^^s — Bd^s mixing and for B decays with the 
underlying quark transitions h suu and h sec. The implication of the 2-3 
symmetry in the up quark sector and on the CKM matrix is studied. This model 
generates FCNC effects in the leptonic sector and predictions for the decay r fill 
and r — >■ iiqq are made. Finally, a deviation of about 2 degrees from the maximal 
value for the atmospheric mixing angle 623 is predicted from the breaking of the 2-3 
symmetry in the charged lepton mass matrix. 

We begin in Sec. 2 with a description of the two Higgs doublet model with 
2-3 interchange symmetry and study the effects of the model in various B decays. 
In Sec. 3 we study the implication of the model for the up quark sector and the 
CKM matrix. In Sec. 4 we study the effects of the model in the lepton sector, 
concentrating on FCNC effects in the r decays and the neutrino mixing matrix. 
Finally in Sec. 5 we conclude with a summary of the results reported in this work. 



Two Higgs doublet models (2HDM) have been studied widely and a particularly in- 
teresting version with FCNC was studied in detail in Ref. JH]- The model presented 
here also has FCNC effects but its origin and structure are different from that in 



We start with the discussion of our model in the the quark sector. Consider a 
Yukawa Lagrangian of the form 



where 0i, for z = 1, 2, are the two scalar doublets of a 2HDM, while Y^j' and S^j' 
are the non-diagonal matrices of the Yukawa couplings. 

For convenience we can choose to express 0i and 02 in a suitable basis such that 
only the Y^'^ couplings generate the fermion masses. In such a basis one can write. 



2 The Model- Quark Sector 



Ref. [IHl. 





(2) 
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The two Higgs doublets in this case are of the form, 

1 f \,^fV2x' 



f V2H^ 
^2 \H^ + iH^ 



(3) 



In principle there can be mixing among the neutral Higgs but here we neglect 
such mixing. We assume the doublet 0i corresponds to the scalar doublet of the SM 
and to the SM Higgs field. In addition, we assume that the second Higgs doublet 
does not couple to the up- type quarks(S'^ = 0). For the down type couplings in 
Eq. HI we have, 

= YgQ,^L<PiD,,B^ + S^Q,,L<p2Dj,R + h.c. (4) 
We assume the following symmetries for the matrices and S^: 

• There is a discrete symmetry under which c?/, ^ — > —d^ ji 

• There is a s — 6 interchange symmetry: s ^ b 

The discrete symmetry involving the down quark is enforced to prevent s ^ d 
transition because of constraints from the kaon system. It also prevents h d 
transitions since Bd mixing as well as the value of sin 2(3 measured in B^{t) J /ipKs 
are consistent with SM predictions. Although there may still be room for NP in 
b ^ d transitions, almost all deviations from the SM have been reported only in 
h s transitions and so we assume no NP in 6 — transitions in this work. 

The above symmetries then give the following structure for the Yukawa matrices. 



Y 



D 



s 



D 



/ yii 








[l 


2/22 


2/23 




2/23 


2/22 


I Sii 











■S22 


■S23 


V 


■§23 


■S22 



(5) 



The down type mass matrix, is now given by = -^Y^. The matrix 

Y^ is symmetric and choosing the elements in Y^ to be real the mass matrix is 
diagonalized by. 



= f/^M^f/=| :^(2/22 - 2/23) 





= -75 75 • (6) 
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It is clear that the matrix U will also diagonalize the matrix when we transform 
the quarks from the gauge to the mass eigenstate via dL,R —>■ Udi^R. Hence there 
are no FCNC effects involving the Higgs 02- 
The down quark masses are given by, 



rrid 



V 



±-^{y22 - 1/23), 
±^(2/22 + 1/23)- 



(7) 



Since << mb there has to be a fine tuned cancellation between 2/22 and 1/23 to 
produce the strange quark mass. Hence, it is more natural to consider the symmetry 
limit ?/22 = 2/23 which leads to = 0. We then introduce the strange quark mass 
as a small breaking of the s — h symmetry and consider the structure. 



D 





2/22(1 + 22;) 
2/22 



with z ~ 2ms/ nib being a small number. Note that we do not break the s — h 
symmetry in the 2 — 3 element so that the matrix remains symmetric. This 
down quark matrix is now diagonalized by. 



ML, 



W 



(±^,yu 
±^^2/22 









±^^(2 + ^22 7 



/I 













(9) 



^V-^2"^ 2' 

In Eq. El we have dropped terms of O(z^). The down masses are now taken to be, 

rris = ±-j=zy22, 

rrih = ±-^(2 + z)y22. (10) 

We find z ~ ±2ms/mi) and now the transformation to the mass eigenstate will 
generate FCNC effects involving 02. For definiteness we will choose the positive 



5 



sign for z though both signs are allowed. The matrix in the mass eigenstate 
basis now has the form, 



S 



D 



s 



D' 





(^22 - S23) 
S23Z 




(11) 



There is now FCNC involving a 6 — * s transition which is proportional to z ^ 
^s/^b ~ where A is the cosine of the Cabibbo angle. 

It is also possible to find a that breaks the s — b symmetry. For example, we 
could choose. 



S 



D 





V 





S22 
S23 





S23(l + 2e) 
S33 



(12) 



with e being a small quantity of the same size or smaller than z. A fit to the 
B Kn data, presented later, will rule out small z/e. The transformation to the 
mass eigenstate leads to the general parametrization. 



S 



D 



where 



Sde 



S 



i4>dd 



D' 



Sde 
















se 



«<Ps6 



se 



+ pe 



She 



pe 

i4'bb 



i^sb 



(13) 



See 



i<t>3 



'11 N /I 

-S22 + -S33 - 523) + [7;S33 



S22)Z - S23e, 



se 



l<Psb 



pe 



/I 1 ^ /I In 

(-S22 + -S33 + S23) + (2-522 - -jS-iz)z + 5236, 

/I 1 N 

(2-^33 - 2-522) + S23Z, 

323^- 



(14) 



There is now an additional FCNC involving b s transitions whose source is the s—b 
symmetry breaking in . Note that the 2-3 off-diagonal elements in contain 
a part that is symmetric under s — b interchange and a part that is antisymmetric 
under the s — b interchange. The parameters in can be obtained or constrained 
from a fit to -B decay data. 

We will not consider this general case in the paper but to make our model 
predictive we will assume a simplified ansatz for in Eq. ^| Here we will consider 
the following structure for the matrix with a small s — b symmetry breaking. 



S 



D 






±s(l + 2e) 




(15) 
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where s and e are real numbers and we have introduced a universal weak phase (j). 
On moving to the mass eigenstate we obtain, 




^ 5^ = e"^ Ts{l + e) ±s{z - e) . (16) 



Here the coupling of the Higgs to the d and the s quark are the same, up to a sign, 
to a good approximation since e is a small parameter. The FCNC 6 — > s transition 
arises from two sources. The first, represented by 2; ~ 2ms /mi,, comes from the 
breaking of the s — b symmetry in the matrix (Eq. |H)) and is symmetric under 
s — b interchange while the second represented by e < 2; , comes from the breaking 
of the s — b symmetry in the matrix Eq. [T3j) and is antisymmetric under s — b 
interchange. 

The size of the matrix elements in are not known. One might argue that 
the size of the matrix elements should be the same as the matrix elements in the Y 
matrix that gives mass to the down quarks. Hence the elements in are given as 

~ y ~ ruq/v where q = d,s,b and v ~ 246 GeV and are small. However it is 
possible that the down quark masses get their masses from a different Higgs boson 
than the top quark jTHI, in which case the vev of 0i giving mass to the down quarks 
can be small, much less than v , thus allowing for larger values of Y and S^. Note 
that it is not necessary for Y and 5*^ to be related and thus we assume no relation 
between them. 

The lagrangian describing the interaction of the Higgs, ifi,2, is given by, 

A = ^JiH, + tH,), 

j"'-' = ^,,(J,(l + 75)^^,±5'*.(J,(l-75)c^^. (17) 
Specifically for b sqq{ q = d,s) transitions the relevant currents, Jij are, 

Jfb'" + Jt" = Ss>,-s{i + i,)b±s;x^-^,)b, 

J?,''' = 5,,g(i + 75)g±5;g(i-75)g. (18) 

Using Eq. (fT?)|) we have, 

Ssti^bs) = ±se'^{zTe), (19) 

Sm = -Sss = ±se"t'. (20) 

After integrating out the heavy Higgs boson. Hi 2, we can generate the following 
effective Hamiltonian with four quark operators. 
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Gf 2ms , i 

H:j, = V,,%'^,^{±0,,T0nn + 0,n-0nL], 

Orr = e*2<^s(l + 75)&g(l + 75)g, 

Oll = e-^20^(l - 75)&g(l - 75)g, 

Orl = s{l + 75)6g(l - 75)g, 

Olr = s(i-75)&g(i + 75)g, (21) 

where rjgg = 1(— 1) for the choices in Eq. (fT6|l . g is the weak couphng and the plus 
(minus) signs in front of Oll,rr correspond to the Higgs exchange Hi 2- Note that 
a new weak phase is associated only with the operators Oll rr for the ansatz in 
Eq. El The Higgs effective potential will be assumed to make one of the neutral 
Higgs lighter than the other and we will be interested in the phenomenology of 
the lightest Higgs only which we will generically denote as H and for the sake of 
simplification and clarity we choose H = Hi in Eq. ()21|). Possible FCNC effects 
associated with the other neutral Higgs boson will be neglected by taking its mass 
to be sufficiently heavy. An effective Hamiltonian will also be generated by the 
exchange of the charged H~^ Higgs boson which will generate transitions of the form 
b cqq. We will not consider such effects in this work. 

Now from the structure of H'^/^ we can calculate nonleptonic decays where hints 
of deviations from the SM have been reported. We will make a semi-quantitative 
analysis of NP effects in the various decays and a more thorough investigation will 
be carried out in a later work ^Tj. We will first fix the parameters of the model 
from the fit results to the Ktt system obtained in Ref. [7]. We will then make 
predictions for several nonleptonic B decays. For the purpose of the paper we will use 
factorization to calculate nonleptonic decays. The impact of non factorizable effects 
will be addressed later in the paper. For the decay B P1P2, where Pi^2 are any 
final state particles, the matrix element of any operator in Eq. (j2ip . in factorization, 
has the structure < O >~< Pi|s7a&|-B >< -P2|^7bQ'|0 > where 7^_b = (1 ± 75) if 
the meson P2 contains the qq component in its flavour wavefunction. An example of 
such a decay is B~ K~tt^ with Pi 2 = K~{tt^). On the other hand for the decay 
B~ K^Ti~ the operators in Eq. |^ have to be given a Fierz transformation. We 
can define the Fierzed operators, 

OIl = -^e-2'^gll-75)&^(l-75)g-^e-2V,.(l-75)&^a'^'^(l-75)g, 
O^RR = -^e%(l+75)te(l + 75)g-^e%-a,.(l+75)fo-^"^(l + 75)g, 
OIr = -^g7M(l-75)fe7'^(l + 75)g, 
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= -^97^(1 + 75)&s-7'^(l-75)g, (22) 

where we have done also a color Fierz and dropped octet operators that do not 
contribute in factorization. 

We can now look at the various nonleptonic B decays and we start with those 
decays with the underlying quark transition b —>■ sdd. 

• B ^ Kit decays: 

Let us denote by A^^ the amplitude for the decay tt^KK In the SM they 

are described to a good approximation by a "tree" amplitude T', a gluonic 
"penguin" amplitude or P', and a colour-favored electroweak penguin (EWP) 
amplitude -P^w In i? — t^K decays, there are four classes of NP operators, 
differing in their color structure: Sa^ihaqpT jqp and SaTibpqjsTjqa {q = u,d). 
The matrix elements of these operators are then combined into single NP 
amplitudes, denoted by A!''^e^^i and ^"^''?e**<' respectively with q = u,d. In 
the presence of NP the most general B ttK amplitudes take the form [71ll6j. 
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v^A°° = -P' - P;^. + ''e'*' + A'^^'^e'"'^ , (23) 
where A''^°"'''e''^' = -^''"e**"+^'''^e**d, A""^"", and A'^^"^ are the NP amplitudes 

m- 

It was found in Ref. [7j that a good fit to the data can be obtained with the 
NP amplitudes \A''''°"'''/T'\ = 1.64, $' ~ 100°, A""'"" ^ 0, and A"'''^ ^ 0. It is 
clear from the the structure of Eq. |^ that A'^'^ = as the NP involves only 
down type quarks. 

Now one obtains, using Eq. |^and factorization, the matrix element relations 
for the decay B~ —>■ K^tt^, 

<0ll{01l)> = -<OnR{O^R)> 

<Ol/?(OL)> = - < OuLiO^r.) > ■ (24) 
These relations follow from the fact that in factorization, 

<Oab> ~ < K'\s-fAb\B >< TT^ld'yBdlO >, 
where 7^^^ = (1 ± 75). We then have, 

<Oll> = <i^-|s(l-75)&|5><7rV(l-75)rf|0>, 
= ~ < K^\sb\B ><n^\d-f5d\0>, 
= - < K-\s{l + 75)6|P >< 7rV(l +75)c?|0 >, 

= - < Orr > . 
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Similar arguments lead to the other relations in Eq. |^ Using the matrix 
element relations in Eq. |21] one obtains, 

< K-n'\H'Jf\B- > = 1^''™™^^*' = ^A,d [2i sin 20x. - 2 cos 20Xa + 2xa] 



Xs 



2ms 

mb m^j ' 
m^ 



w_ 

Xa C 9 2 ) 

Am = < K-\sh\B- >< 7r°|(i75(i|0 > . (25) 

To make estimates we choose e ~ z = 2ms /mi, and cos 20 ~ — sin 20. Choos- 
ing 20 = —45° and the fact that tan$' ^ cot0 leads to $' ~ 113°, which is 
consistent with the fit obtained in Ref. 0- We observe that if Xa/Xs = e/-2 ~ 
then $' ~ 90° which is also consistent with the fit obtained in Ref. ,7J. How- 
ever, with Xs/Xa = z/e ^ we obtain $' ~ which is inconsistent with the 
fit in Ref. [7j. Hence only e < 2 is allowed as indicated earlier. Turning now 
to the amplitude A"^''^ we have, 

<K\~\H'f^B-> = l^'^-'^le 

1 Gf 

[2Frfd(zsin20Xs - cos20Xa) + ^XaGdd] , 



2N,V2 

Fdd = <n-\db\B- ><K'>\s-f5d\0>, 

Gdd = <n-\d^^b\B- ><K^\sYl5d\0> . (26) 

It is clear from Eq. that A'^''^ is suppressed relative to ^"^'^omfc j2W^^ 
and hence small, which is again consistent with the fit obtained in Ref. [7]. 
Now using \A'''^°'^^/T'\ = 1.64 we find, using naive factorizationjlSj . 

_ 2fc sin 2(f)XsAdd 

V:bVus (7r°| 67/.(l - 75)m \B-) {K~\ u^il - l,)s |0) 

k = v^l + tan2 0- (27) 

One can then convert this to 

■ U{ml-ml)F^/V2 / ^C2\,,.,,. , 

2k\ sm20 Yc = 7= 1-64 ci H K.^Ks 

' "^'-"^ [{ml - mi,)/(m, - ms)]F,^iml/2md)U/V2 V ' iV J ' 

(28) 

We take {U/ U)iF^ / F^^) ~ 1, \V:,Vus/V,lVts\ = 1/48 and Ci + c^/N, = 1.018. 
Taking the masses from the Particle Data Group [2111, we find, 

Xs = TT^^i^-O^T^lVtlVtsl (29) 
k\sm2(p\ bMeV 
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This then leads to, 



m 



w 



nib 



2s^A;| sin 201 



O.OSSlV^lVtsl 

o-omiVts\ 



md = A MeV, 
md = 8 MeV, 



(30) 



and finally to. 



mw 



V2s 



9 



rris fc|sin2</)| \/2s rUs fc|sin20| 



\/2s 



k \ sin 20 1 



A.lMw(- 
9 



< niH <6.SMw(-), 
9 



(31) 



where we have used = 100 MeV and nib = 5 GeV. For s ~ 1 we obtain 
niH ~ 600 - 900 GeV. 

Our estimate of the parameters of the model, like the new physics weak phase 
and the higgs mass mn, used the factorization assumption to calculate the 
nonleptonic amplitudes. It is therefore important to address the impact of 
non factorizable effects on the values of the parameters of the model. We note 
that to extract the parameters of the model we have to calculate the ratio of 
the tree amplitude in the standard model relative to a new physics amplitude. 
Since non factorizable effects in the tree and the new physics amplitudes are 
expected to be small [Tallin] we do not expect these effects to significantly alter 
the values of the parameters of the model extracted using factorization. A more 
precise determination of the model parameters using QCD factorization 
will be carried out in a future work [T7j . 

B K*p decays: 

Here the final states contain two vector mesons and so from angular momentum 
considerations there are three amplitudes. They are usually chosen to be 
the longitudinal amplitude, Aq, and the two transverse amplitudes Ay and 
A±. In the SM, Ao is the dominant amplitude and the transverse amplitudes 
are ©(l/m^) because of the V — A structure of the weak interactions. Now 
it is clear from the structure of the operators in Eq. (PT|) that the scalar 
operators do not contribute to vector- vector final states in factorization. Hence 



K- 



On the other hand 



K^*p gets contribution from the the Fierzed operators in 



we expect no new physics effect in the decay B 
the decay B~ 

Eq. 1221 The tensor operators in Of^ contribute to the transverse amplitude 
in the leading order in the large limit unlike in the SM. This follows 
from the fact that the matrix elements of the tensor operators contain the 
piece ~< K*'^{q^)\safj,^d\0 >~ {e^qu — e^q^j) that contribute to the transverse 
amplitudes in the large tub limit. For a longitudinally polarized state the 
polarization vector can be approximated as, ~ neglecting (7^)^ effects, 
and hence the contribution to the longitudinal amplitude, Aq, arise only at 
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0{mv /TnB). The operators Of^B.RL contribute in the leading order in the large 
iTiB limit to the longitudinal amplitude, Aq while the transverse amplitudes 
are suppressed by O(^) PCSl. 

Hence our prediction is that the decay K^*p^ should be dominantly 

longitudinally polarized as in the SM because there is no new physics effect 
in this decay, while the decay B^ ^ K^*p~ may have a sizeable transverse 
polarization. These predictions appear to be consistent with present exper- 
iments where BaBar measures the longitudinal polarization fraction, fi, for 
B- K-*p^ to be O.gelHs 

± 0.04 while the BaBar and Belle measurements 
of /l for B- K^*p- are 0.79 ± 0.08 ± 0.04 and 0.43 ± 0.11°-[]5o2 respec- 
tively giving an average Jl of 0.66 ± 0.07, thereby showing sizeable transverse 
polarization in this mode jl2[ [13\. 

To calculate specific quantities we again use factorization to compute the non 
leptonic amplitudes. We do not expect non factorizable effects to significantly 
alter the pattern of new physics effects in the B decays as our conclusions to 
a large extent rely on heavy quark theory considerations. As far as specific 
values of various quantities, like polarization fraction, CP violation etc, in the 
B decays are concerned it is possible that non factorizable effects will change 
the values of these quantities. Hence, our calculation of these quantities should 
be taken as rough estimates. A more precise estimate of the quantities will be 
carried out in a future work|17j. 

For the A± amplitude the factorized matrix element satisfy < Of^^^ >=< 
Orr > while for the Ay amplitude the factorized matrix element satisfy < 
Of^ >= — < O^j^ > and so one obtains, 

A± ~ [Xs cos 20 - ixa sin 20] , 
^11 ~ [~Xa COS 20 + ixs sin 20] . (32) 

Note that even though the Fierzed tensor operators are colour suppressed 
their effects for the transverse amplitudes can nonetheless be significant as the 
SM penguin amplitude for VV(V=vector) final state is smaller than the PP 
(P=pseudoscalar) final states yLB|. It follows from Eq. (jH^ . for Xs ^ Xa, that 
the transverse polarization fractions f± ~ f\\ where /t ~ I^Tp (T =±, ||). We 
also expect no triple product asymmetry from the interference of A± and 
A\\ because the relative weak phase difference is zero (or vr). A triple product 
asymmetry from the interference of the longitudinal amplitude Aq and A± can 
be significant. 

We now consider nonleptonic B decays involving the underlying b sss transi- 
tions. 

• i? — > (pKs, B Tj'Ks decays: We begin with the decay B (pKg where only 
the Fierzed operators in Eq.|221 contribute as is a vector meson. In factoriza- 



12 



tion, for this mode, we have the matrix elements relations ,< Oll >=< Orr > 
and < Olr >=< Orl >■ We can therefore write, 



ASM , aNP , aNP 

^ ^ ^LR±RL ^ ^LL±RR-> 



-^Vt,V*Z 



t t t ^ t ^ t ^ t 

% + 04 + — -a-j — -ttg — 2*^10 



A 



NP 

LR±RL 



aNP 

^LL±RR 



Xs 

z 



1 Gp 



Xa 



Z [xs COS 20 - ixa sin 20] 



rO.OSlKK, 



A;|sin20r'"^' 
2f^m^FBK{ml)e* ■ ps, 



(33) 



where the SM contribution can be found in Ref. |2T] and we have chosen 
md = 6 MeV. 

We can now calculate sin(2/?)e// from. 



sin(2/?)e// 
A/ 



2Im[X 



:i + iA,pi 



-2if3 



A 
A' 



(34) 



where A = A'^^i^ + and A is the amplitude for the CP-conjugate process. 
Note that from Eq. (|3l]p . sin(2/3)ejj is independent of Z and hence free from 
uncertainties in the form factor and decay constants. We also observe that the 
deviation from the SM expectation of sin(2/3)e// comes from which is 

however colour and ^ suppressed, leading to a rough estimate, sin(2/5)e// ^ 
0.43 which agrees well with the present experimental average of 0.44 ± 0.27 ± 
0.05 jHHHEl • Note that the fact A^^^^^ does not carry any new weak phase 
is a consequence of the form in Eq. and in the general case ( Eq. IT^ there 
will be a new physics phase associated with this operator. Furthermore, if 
Xa/Xs = ^/z ~ then there is no new weak phase in this decay and the 
prediction of sin(2/5)e// should be given by the SM. 

The decay B Tj'Kg is more complicated as both the h sdd and h sss 
transitions contribute ^^j- The final states are pseudoscalars and so both the 
operators in Eq. Eland the Fierzed operators in Eq. [221 for b — >■ sss transition 
contribute. Hence in this model this decay could show a significant deviation 
from the SM prediction for sin(2/?)e//. For the "standard" rj' wave function 
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^ ^ {uu+dd_ 2ss) ^j^g j^-p ^ _^ ^ _^ amplitudes add or partially 

cancel depending on the choice in Eq. CSl In the t] they would cancel or 
add producing a much smaller or bigger deviation from the SM in 5 — >^ rjKg 
compared to 5 — rj'Kg. It should be noted that the decay B — > rjKs has a 
much smaller branching ratio, probably because of destructive interference in 
the SM amplitudes [231 • 

• B'^ —>■ (j)K* decays: Here the final states contain two vector mesons like the 
B~ — > p~K* decays and the predictions here are very similar. This decay 
gets contributions only from the Fierzed operators in Eq. |221 The tensor 
operators in O^^ ^p, can contribute significantly to the transverse amplitude 
in the large ms limit [13] thus explaining the low longitudinal polarization 
fraction, fi = 0.45 ± 0.05 ± 0.02 [TTl Ej, measured in this decay. Given the 
large NP weak phase we could also observe a sizeable triple product asymmetry 
and/or direct CP asymmetry in this decay. 

To summarize, the model can provide explanations for the deviation from SM 
seen in several rare B decays. It also makes specific predictions in other B decays. 
For example, decays going through b suu transitions like Bg K^K~ should not 
be affected. The decays going through h sec transitions like i?^ D^^*^D^^f*^ 

should also not be affected and so this decay along with Bd — > can be 

used to measure the angle 7 without NP pollution. The model also has no effect 
in the decays B^^s J/V^-ft''-*-* (0, r]'). A detail study of the predictions of this model 
in B decays will be presented elsewhere ^17) . 

We now discuss the important case of B^^s ~ Bd^s mixing. By construction of 
the model there is no effect in Bd — Bd mixing. In B^ — Bg mixing this model 
can produce a new CP phase through contributions associated with the operators 
Bll = b{l — 75)^6(1 — 75)5 and Bjm = 6(1 + 75)56(1 + 75)5. In the vacuum inser- 
tion approximation, the matrix element of Bll and -B_rr are the same as only the 
pseudoscalar currents contribute. The contribution to Bg mixing comes from the 
sum < S*le~'^'^'^ Bll + Sl^e^'^'^ B rli > which contain a term ~ iXsXa sin 20 which is a 
source of new weak phase in Bg — Bg mixing. This term will be small or vanish if 
Xa/Xs = e/ z ^ 0. Hence the presence of a new weak phase in the Bg — Bg mixing 
requires the breaking of the s — b symmetry in the Yukawa coupling of the second 
Higgs doublet. 

Finally, we note that the model will produce new effects in 6 — sl^l~ and 
Bg decays which will depend on the couplings s and si of the second Higgs 

doublet to the quarks and leptons. A detailed study of such processes in the model 
will be discussed elsewhere [T7]. 
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3 Up Quark Sector 



So far we have neglected the up sector and it is possible that there will be FCNC 
decays in that sector also. However the diagonalizing of the up quark mass matrix 
is connected to the CKM matrix via Vckm = Vl'^ where W in Eq. El diagonalizes 
the down quark matrix and transform the left handed up type quarks from the 
gauge to the mass basis. If we assume the up sector has the same symmetry as the 
down sector then Vl ~ W. We can write, 



^(1 + I^J 

- 2^«) 


72(1 - 1^'^) 





(I 







= 





-71(1- 


- 












/I 












-75(1- 


- \zd 






71(1+ 


\zd) 



(35) 



where Zu ~ ±2mc/mt and Zd ~ ±2ms/mb break the 2-3 symmetry in the up and 
the down sector. The CKM matrix is now obtained as, 



CKM 





1 _|_ Zd 

) ( ^tt 







Zd) 



\ 



(36) 



We see that we get the right order and sign for the CKM element Vd, and Vts. To 
obtain the realistic Vckm we have to introduce the Cabibbo angle A as a symmetry 
breaking effect. Possible FCNC effects in the top sector will be discussed in a later 
work[T7]. 



4 Leptonic Sector 

In this section we study the consequences of the 2 — 3 symmetry applied to the 
charged lepton sector. We assume that the structure for the charged lepton mass 
matrix, is given by, 

(hi \ 

Y"^ = \ Q l22{l + '2zi) I22 , (37) 
V k2 I22 1 

where zi ~ 2m^/m^ is a small number. This charged lepton mass matrix is now 
diagonalized by. 











±^,z,l 









Mi,^ = Wi^M^'Wi = I ±^,zih2 



\ 


±^(2 + ^0^22 7 
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w, 



/I 









h^"^ 2 



Just like the down quark sector we have broken the ^ 
In the symmetry hmit, 



/I 



W, 



U, 









1 


1 


V2 


V2 


1 


1 


V2 


V2 



(3^ 



r symmetry by the /i mass. 



(39) 



The Yukawa interaction associated with the second Higgs doublet is now taken 
similar to Eq. as, 




(40) 



where si and ei are real numbers and we have introduced a universal weak phase (pi. 
On moving to the mass eigenstate we obtain, 











^Si(l + eO ±si{zi-ei) 
\0 ±si{zi + ei) ±si{l + ei) 



(41) 



To simplify our discussion we assume e/ ~ 0. The model now generates FCNC 
interactions of the r like r — > nil where / are muon or electrons. Note that decays 
r — s> filJb where la 7^ h are forbidden. The effective Hamiltonian for such decays is 
easily written down as, 



Oll 

Orr 

Olr 

Orl 



1 



q2 



4 m-j- m jj 



[±Oll ± Orr + Olr + Orl] 



e-^^*'/i(l-75)r/(l-75)/, 
e^2<^'/i(l + 75)^1(1 + 75)^ 
/i(l - 75)rr(l + 75)/, 
/i(l + 75)rr(l - 75)Z. 



(42) 



To further simplify the discussion we will choose 0; = 0. In ref. [2S] the FCNC 
leptonic transitions were studied in an effective operator formalism where the coeffi- 
cient of the four quark operators were taken to be ~ 47r. This lead to a constraint on 
the scale of NP, A ~ 10 TeV. To compare to our model we have the correspondence. 



47r 

A2 



4 rrir rri^j 



(43) 
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Choosing ~ 1 we find niH ~ 340 GeV. From B decays we found thh ~ 600 — 900 
GeV and so we predict tlie brancliing ratio of r fill to be below the experimental 
bounds [20] by about (340/900)^ - (340/600)^ ~ 2 x 10-^ - 10"^ A careful analysis 
would include varying the parameters (pi and si and will be presented elsewhere [T7]. 
We can also study the decay r fiqq which can be studied in the decays r — > yU7r*^(p). 
The experimental bounds produce similar size limit on the Higgs mass as the decays 
T — > fill. Note that for 0; = only H2, which is a pseudoscalar, can contribute to 



r 



fiTT and the contribution from Hi exchange vanishes. 



Let us now turn to the neutrino sector. The neutrino mixing, the PMNS matrix, 
arise from the lepton mass Lagrangian as follows: 

Cm = vlC~^My,a^pu + e^^LMlpeR + h.c. (44) 

Diagonalizing the mass matrices by the transformations U]^ M-yUy = J^^iag and 
UjM'^V = MJj^^, one defines the neutrino mixing matrix as Upmns = UjUu- We 
will parameterize Upmns as follows: 



U 



PMNS 



( C12C13 

-S12C23 - Ci2S23Si3e* 



S12C13 
C12C23 - Si2S23Sl3e*' 



V ■512S23 - Ci2C23Sl3e*'^ -C12S23 - Ci2C23Si3e*' 




(45) 



where K 



diag{l, e 



The 2-3 symmetry in the leptonic sector leads to a simplified form of the PMNS 
matrix, with S13 = 0, given by. 



/ 



PMNS 



C12 

1 




(46) 

In the basis where the mass matrix of the charged leptons is diagonal, the left- 
handed flavor-based (symmetric) neutrino mass matrix m! is related to its diagonal 
form, My^diag = diag[mi,m2,ms] by. 



u 



PMNS^^-^v^PMNS 



The form of Up^NS J^^^ follows from the fi — r symmetry in Ml p]. 
Note that we can express Up^j^g, using Eq. [HHl as. 



u 



PMNS 



(47) 



(4J 



where 



Ul 









1 






75 


1 








-Si 




Cl 


;) 








(49) 
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So the neutrino matrix, is just a combination of a simple rotation matrix and 
a phase matrix. Now the breaking of the 2-3 symmetry in the leptonic sector will 
change Uj and Upmns to, 



U 



w. 



/I 



vo 



PMNS 



Cl2 




V2 



V2^ 



(50) 



V 7l(l + tK^ -775^12(1 + 1) ;l5(l-f 



where zi = ±2^. This corresponds to S23 = :^(1 + f") C23 = 775(1 — f-) and to 



V2^ 



an atmospheric mixing angle of 623 ~ 43.26 for the negative sign of zi. 

Finally, we point out that this model will have interesting collider signatures. 
The heavy Higgs, H, according to the structure in Eq. UHl couples equally to all 
three generations, to a very good approximation, unlike the usual SM Higgs which 
has couplings proportional to the mass. The effect of electroweak precision mea- 
surements and collider signatures will be explored elsewhere jT7]. 



5 Conclusions 

In summary, there are now several B decay modes in which there appear to be 
deviations from the SM predictions. These deviations could signal the presence of 
beyond the SM physics. In this work we were interested in a NP scenario that can 
provide a single solution to all the deviations. We considered a two Higgs doublet 
model with a 2-3 symmetry in the down type quark and the charged lepton sector. 
The breaking of the 2-3 symmetry, introduced by the strange quark mass and the 
muon mass lead to FCNC in the quark sector and the charged lepton sector that 
are suppressed by ^ and ^ in addition to the mass of the heavy Higgs boson. 
Additional FCNC effects of similar size were generated from the breaking of the 
s — b symmetry in the Yukawa coupling of the second Higgs doublet. From a fit 
to the B Kit data we found the mass of the lightest neutral Higgs in the second 
doublet to be of the order niH ~ 600 — 900 GeV. We made several predictions in B 
decays listed below: 

• A sizeable transverse polarization in the decays B^ —>■ (j)K* and B~ —>■ p~K^* 
was predicted but not in B~ p^K*~ . This is consistent with present 
measurements jT2| IH] • We also predicted the possibility of observing a size- 
able triple product asymmetry and/or direct CP asymmetry in the decays 
5^ ^ (l)K* and B' p-R"^*. 

• The sin 2/3 measurement in B^{t) (pKs was predicted to show a small devi- 
ation from the SM value( if Eq. [121 is assumed) but it was found that sin 2/3 
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measured in B^{t) rj'Ks could have significant deviation from tlie SM pre- 
diction. 

• Tlie decays witli tlie quark transition h — > suu and h — ^ sec were found to be 
unaffected 

• We found a new source of weak pliase from our model in Bs — Bg mixing while 
Bd — Bd mixing was not affected. 

We studied the implication of the 2-3 symmetry extended to the up sector and found 
that we could generate the right sign and size of the CKM matrix element Vch and 
Vts- We then studied FCNC effects in the lepton sector. The lepton flavour violating 
decays r — ^ were found to be below the present experimental limits by a 

factor of 2 X 10^^ — 10~^. Finally, we found that the breaking of 2-3 symmetry in 
the lepton sector could lead to deviations of the atmospheric neutrino mixing angle 
from the maximal value by ~ 2 degrees. 
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